VECTOR INTERPRETATION OF THE MATRIX 
ORTHOGONALITY ON THE REAL LINE 



A. BRANQUINHO, F. MARCELLAN, AND A. MENDES 



Abstract. In this paper we study sequences of vector orthogonal polynomials. 
The vector orthogonality presented here provides a rcintcrpretation of what 
is known in the literature as matrix orthogonality. These systems of orthogo- 
nal polynomials satisfy three-term recurrence relations with matrix coefficients 
that do not obey to any type of symmetry. In this sense the vectorial reinter- 
pretation allows us to study a non-symmetric case of the matrix orthogonality. 
We also prove that our systems of polynomials are indeed orthonormal with 
respect to a complex measure of orthogonality. Approximation problems of 
Hcrmite-Pade type are also discussed. Finally, a Markov's type theorem is 
presented. 



1. Introduction 

In the late eighties of the last century, the following problem attracted the 
interest of many researchers. 

When a sequence of monic polynomials, {p n }n£N, satisfying a recurrence relation 

N 

X N p n (x) = C nfi p n {x) + ^ [c n ,kPn-k(x) + C n+k . k p n+k (x)] , (1) 
fc=l 

where c n fi (n = 0, 1, . . .) is a real sequence and c ntk , (n = 1,2, . . .) are sequences of 
complex numbers for k — 1, 2, . . . , N with c^jv 7^ 0, is related with some kind of 
orthogonality? 

Several authors (A. J. Duran, F. Marcellan, W. Van Assche, and S. M. Zagorod- 
nyuk, among others) were interested on this subject. Their contributions re- 
vealed an enormous interdisciplinarity between different kinds of orthogonality (like 
Sobolev orthogonality, orthogonality on rays of the complex plane) and several ap- 
plications, mainly quadrature formulas. From the extensive bibliography on the 
subject we stand out the references PQ [31 HI [11 [HI [13 [13 [III [IS [T51 [T71 US] - 

In his work [3], A. J. Duran presents for the first time a Favard's theorem 
for sequences of polynomials {p n }neN satisfying recurrence relations like Few 
years later this result was reformulated by the author together with W. Van Assche 
in [6] where they stated the connection between sequences of matrix orthogonal 
polynomials and sequences of polynomials that satisfy a higher order recurrence 
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relation. As an application, the authors gave an interpretation of a Sobolev discrete 
inner product. 

In this context the authors considered a positive integer number ./V and the op- 
erators i?Ar, m , m — 0, 1, . . . , N ~ 1, defined on the linear space of polynomials, P, by 

i.e, the operator RN,m takes from p just those powers with remainder m (mod- 
ulus N) and then removes x m and changes x N to x proving then, the following 
result. 

Theorem 1 ( L 3J). Suppose that {p n }nen, with degp n = n, is a sequence of poly- 
nomials satisfying a (2N + 1) — term recurrence relation as ([T|) and let {P n }neN be 
a matrix polynomial sequence defined by 

RN,o{PnN){x) ■■■ RN,N-l(PnN){x) 

P n (x) = : ••. : 

_RN,o(P(n+l)N-l)(x) ■■■ RN,N-l(P(n+l)N-l)(%)_ 

Then, this sequence is orthonormal on the real line with respect to a positive defi- 
nite matrix of measures and satisfies a three-term recurrence relation with matrix 
coefficients. 

Conversely, suppose that {P n }nefh with P n = (P™^)m~j=o> * s a sequence of or- 
thonormal matrix polynomials or, equivalently, they satisfy a symmetric three-term 
recurrence relation with matrix coefficients. Then the scalar polynomials defined by 
N-l 

PnN+M = J2 x j Pn' j {x N ) (« e N, < m < N - 1) 

3=0 

satisfy a (2N + I) — term recurrence relation of the form |T]). 

Taking into account the current relevance of the subject our work is concerned 
with the analysis of higher order recurrence relations, in this case, of order 2N + 1 

2JV-1 

h(x)p n (x) = cf^-'pn+Nix) + Y, <XN-\-KPn + N-X-klx) (2) 

fe=0 

where h is a polynomial of fixed degree N and where c" +Ar_1 , n > 0, are real 
sequences for j = n — N, . . . , n + N — 1 with c™^ -1 an d initial conditions 
on pi for i = 0, . . . , N — 1 are given. 

We begin by pointing out that in the structure of the recurrence relation ^ the 
polynomial h is a generic polynomial with fixed degree N and their coefficients do 
not satisfy any kind of symmetry. 

Our aim is to analyze this more general case by studying the sequences of poly- 
nomials satisfying such a kind of recurrence relations in order to find out what type 
of orthogonality is associated with them. On the other hand, as an application, we 
expect to obtain some known results. 

Let us consider the family of vector polynomials = {[p\ ■ ■ ■ pn] T ■ Pj eP}, 
and Mat x jv(M) the set of N x TV matrices with real entries. Given a polynomial h, 
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with dcg h = N, we can split the linear space of polynomials, P, using the basis 

{1, x, ... , x N -\h(x),xh(x), x N - 1 h(x),h 2 (x),xh 2 (x), . . .}. (3) 

Then, let {TjljeN be a sequence of vector polynomials such that fj(x) = (h(x)) 3 
7 (x), where %(x) = [lx ■■■ x N ^] , j G N. Let {p m } m£ N be a sequence of 
polynomials, degp m = m, m S N. We define the associated vector polynomial 
sequence {T> m } meN by 

■Bm = [p m N ■ ■ ■ P{m+l)N-l] T , n E N . 

A scalar polynomial p m N+k of degree mN + k, with < k < iV — 1, can be 
expanded in the basis ((3|) as follows 

m JV-1 

p m AT +fe (a;) = ^ dijx j h l (x). 

If we consider the operator Rh,N,j that takes from p m N+k the terms of the form 
a^jX 3 h l {x) and then removes the common factor x J and change h(x) to x, we get 

N-l 

PmN+k(x) = ^ X 3 Rh,N,j(PrnN+k)(M x ))- 
j=0 

It is easy to see that we can write !B m in the matrix form 

•B m {x) = V m (h{x))T (x), (4) 
where V m is a N x N matrix polynomial of degree m given by 

Rh,N,o(PnN)(h(x)) ■■■ Rh,N,N-l(PnN)(h(x)) 



V m {h{x)) 



Rh,Nfi{P{n+l)N-l){h(x)) ■■■ Rh,N,N-l(P(n+l)N-l){H x )). 



and ?o(x) = [l x ■ ■ ■ x 1 . Equivalently we can write the elements of the se- 
quence of matrix polynomials {VmjmeN in the form 

m 

v m (h(x)) = J2 B T( h ( x )) j > 

j=o 

where (B" 1 ) is a family of matrices with real entries. 

First we want to prove that if a sequence of scalar polynomials {p n }neN satisfies 
a recurrence relation like @ then there exists a sequence of vector polynomials 
denoted by {® m }„ l£ N and a sequence of matrix polynomials {Kn}meN defined by ^ 
that satisfies a recurrence relation with matrix coefficients and the converse is also 
true. 



Notice that we can rewrite ^ changing n by n + N — 1 , 

2JV-1 

, / s / s n+2(N-l) , \ »+2(iV-l) / n 

h{x)p n+N - 1 {x) = C„ +2 ] v _ 1 'p n+2 N-l(x) + C n+ 2{N-V)-kPn+2{N-l)-k{x) (5) 

fc=0 

and then, consider the N equations associated with © and ([5]). 
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A straightforward calculation yields that the above system of N linear equations 
can be written in the matrix form 



h(x) 



Pn{x) 
Pn+N-l{x)_ 

Pu-n{x) 

_p n -l(x)_ 

Introducing the change of index n = mN in the above relation we get 

h(x)-B m (x) = A m T> m+1 {x) + B m -B m (x) + C TO 2 ro _i(a;), m > 1, (6) 

Similarly, if we take into consideration Theorem [T] given by A. J. Duran, and 
instead of using the canonical basis for the linear space of polynomials, P, we deal 
with the basis then we get a sequence of polynomials {Kn}meN that satisfies 
the three-term recurrence relation 

V m (z) = A m V m +i{z) + B m V m (z) + C m V m -i(z), m > 1 

with some given initial conditions. 

Notice that if we multiply this last relation by CPo then we obtain the recurrence 
relation for the sequence of vector polynomials {H m }meN given by ([6]). Finally, 
from ([6|) and taking into account the structure of the recurrence relation as well as 
the expression of the vector 53 m , we get the (2iV+ 1)— term recurrence relation (J2J. 

Theorem 2. Let {p n }n£N be a sequence of scalar polynomials, {S TO } TO £N the se- 
quence of vector polynomials with 

BmO) = [PrnN(.%)PmN+\{x) ■ ■ ■ P(m+1)N-1 {%)] T , 

and {V m } m £iq the sequence of matrix polynomials given in (J4J. Then, the following 
statements are equivalent: 

(a) The sequence of scalar polynomials {p n } n eN satisfies ([2|). 

(b) The sequence of vector polynomials {23 m }meN satisfies 

h(x)"B m (x) = A m S m+ i(a;) + B m 'B m (x) + C m 'B m -i(x), m > 1, 

with initial conditions < B_i(x) = Oat x i and r Bo(x) given. 

(c) The sequence of matrix polynomials {V m } me n satisfies 

zV m (z) = A m V m+ i(z) + B m V m {z) + C m V m -i(z), m > 1, 
with initial conditions V-i(z) — Onxn and Vq(z) a fixed matrix. 



Pn(x) 



Pn+N-l(x) 



n-f-JV —l 
n+N 



n+2N-2 
n+N 



n n+2N-2 
-■n+2N—l. 



Pn+N(x) 



Pn+2N-l(x) 



„n-NV-l 



„n+2N-2 



-n+JV-1 



n+2N-2 



'n-JV 



n+n-L 



n+2N-2 
-ii- 1 
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The matrices A m , B m , andC m in the recurrence relations are given, respectively, by 



Jm+1)N-1 
Z (m+1)N 

Jm+2)N-2 
"(m+l)iV 



Jm+2)N-2 
'(m+2)jV-l. 



Jm+1)N-1 
2 mN 



Jm+2)N-2 
2 mN 



and 



(m+i)JV-r 

C (m+l)Af-l 

(m+2)N-2 
C (m+1)N-1. 

(m+l)iV-l 
(m-l)iV 





Jm+1)N-1' 
"mN—1 



Jm+2)N-2 
"mN—1 



Now we consider the sequence of matrix polynomials {V m } m< =n defined by 

xV m (x) = A m V m+ i(x) + B m V m (x) + C m V m -x{x) m > 0, (7) 
with initial conditions V—i(x) = Onxn and Vq(x) — Inxn ■ 

The first question is to know when a sequence of matrix polynomials defined 
by ([7]) is related to the matrix orthogonality. 

If C m — A m _i and B m is a positive definite matrix of measures W supported 
on the real line then the polynomials {V m } m ^ are orthonormal with respect to a 
left inner product, i.e., 



(Vi,Vi, 



V l (x)dW(x)V 1 T (x)=S lJ I NxN . 



(8) 



In the last years several authors have studied analytic properties of matrix or- 
thonormal polynomials (see for example HI El E] ) and their connections with the 
spectral theory of linear differential operators with matrix polynomials as coeffi- 
cients. 

In the case when neither C m = A^ n _ l nor B m are symmetric we cannot guar- 
antee that the system of matrix polynomials {V m } m ^ satisfying the recurrence 
relation {3 is orthogonal with respect to a inner product induced by a positive 
definite matrix of measures W. 

In [2] the authors presented a result that characterizes the existence of a matrix 
of measures W such that the system of polynomials {V m } m ^ is orthogonal in the 
sense of ©■ In fact, if the matrices A m and C m , for m G N, in the recurrence 
relation ([7]) , are non-singular then there exists a matrix of measures on the real line 
with a positive definite Hankel matrix as moment matrix such that the system of 
polynomials {V m } me ^ defined by ([7]) is orthogonal with respect to the measure W 
in the sense of ([5J) if and only if there exists a sequence of non-singular matrices 
{-Rmjmei such that the following relations hold: 

• R m B m R^ is symmetric, Vm € No , 

• R T m R m = C~ T • • • Ci T (rTr q )A ■ ■ ■ A m - X , Vm e N . 

In this contribution we prove that a recurrence relation (J7J characterizes a dif- 
ferent kind of orthogonality. The structure of the paper is as follows: In section 2, 
we present the algebraic theory of the sequences of vector polynomials. In this con- 
text, we define a vector linear functional and we introduce the concept of right and 
left-orthogonality with respect to this linear functional. In section 3, we present 
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a reinterpretation of the matrix orthogonality in terms of the vector orthogonal- 
ity showing that there are two sequences of matrix orthogonal polynomials with 
respect to a matrix of measures, not necessarily positive definite, which are bi- 
orthogonal with respect to a vector linear functional. In section 4, we analyze two 
type Hermite-Pade approximation problems and, finally, a Markov's type Theorem 
is deduced. 

2. Vector orthogonality 

Let (P^)* be the linear space of vector linear functionals defined on the linear 
space of vector polynomials with complex coefficients P N , i.e., (P N )* is the dual 
space of P N . In this space we define a vector of functionals as follows. 

Definition 1. Let u 3 : P — > R with j = 1, . . . , N be linear functionals. We define 
the vector of functionals U = [u 1 ■ ■ ■ u N ] in P N with values in Mat x w(IR), by 

••• (u N ,pi)~ 

'. ' • '. ' 

iy,p N ) ■■■ (u N ,p N )_ 
where "(.)" means the symbolic product of It and T T . 

Let A (x) = Ak x k , where Ak G M7v X Ar(M), be a matrix polynomial and U be 

a vector of linear functionals. Let us consider the vector of linear functionals, the 
so called left multiplication ofli by A, that we will denote by All, such that 

(AU)(y) := (AU.y T f = ^(x k U)(T) (A k f . 

k=0 

We will introduce the concept of sequence of vector polynomials left-orthogonal 
with respect to the vector of linear functionals II and we will prove that It is 
quasi- definite, i.c, there exists a unique sequence of vector polynomials, up to the 
multiplication on the left by a non-singular matrix, that is left-orthogonal with 
respect to It. 

Definition 2. Let {p n }neN be a sequence of scalar polynomials with degp„ = n , 
n E N. Let ft. be a polynomial of fixed degree N, {CB m } me N be a sequence of 
vector polynomials with ® m (x) = [p mN (x) p mN+1 (x) ■ ■ ■ P( m +i)N-i(%j\ , and let 
II = [u 1 ■ ■ ■ u N ] be a vector of linear functionals. {23 m } mG N is said to be left- 
orthogonal with respect to the vector of linear functionals U if 

(a) (h k U) (S ro ) - NxN , k = 0, 1, . . . , m - 1. 

(b) (h m U) (!B m ) = A m , to £ N, where A m is a non-singular upper triangular 
matrix. 

We introduce the notion of moment associated with the vector of linear func- 
tionals U. Taking into account that with 7j{x) = (h(x)y? (x) and 
J"o(a;) = [Ix ■ ■ ■ x N ^ 1 ] T , is a basis in the linear space of vector polynomials P^, 



lt(?) := (U.? T ) T = 
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we denote (x k U)('J'j) = U k the j-th moment associated with the vector of linear 
functionals x k li. 

The Hankel matrices associated with IX are the matrices 

IXo • • • IX m 



D„ 



11 



Uo 



m S N, 



where lij are j-th moments associated with the vector of linear functionals IX. IX 
is said to be quasi- definite if all leading principal submatrices of D m , m 6 N, are 
non-singular. 

The following result provides a necessary and sufficient condition for the exis- 
tence of a sequence of vector polynomials which are left-orthogonal with respect to 
the vector of linear functionals IX. 

Theorem 3. Let IX be a vector of linear functionals. Then IX is quasi- definite if, 
and only if, there exists a unique sequence of vector polynomials { r S> m }meN such 
that 15m = ^2j—Q otf^j, where aj 1 € MjvxaK^) with a™ is non-singular lower tri- 
angular matrix and a unique sequence, {A m } mg pj, of non- singular upper triangular 
matrices such that (h k U) ( r B m ) = A m dk,m, k = 0, 1, . . . , m, m £ N . 



Proof. To prove that IX is quasi-definite. Let {'B m }meN be a sequence of vector 
polynomials with H m = X)j=o a J l -^i' where a" 1 S Mat x at(R) and {Jjj-jeN is a 
basis in ¥ N , such that V^x) = (h(x)) j 'P (x) ) %(x) = [l x ■ ■ ■ x N ' 1 ] T . 

From the orthogonality conditions, the vector sequence of polynomials { < B m } m eN 
is left-orthogonal with respect to the vector of linear functionals IX if, for k = 
0, . . . , m — 1, 



(h k U) (£„-, 

and for all m £ N, 
(h m U) (B r 



(h k U) (EafTj 

3=0 



3=0 



(h m U) ajy 3 ) = £ < (h m U) (IP,) = A, 



3=0 



3=0 



Taking into account (ft. fc lX)(J' : ,) = IX^+fc), the above conditions can be read as 



[ < 



U(3> ) 



U(3> m ) 



= [ 



(9) 



U(9 m ) ■ ■ ■ U(T 2m ) 

For m = 0, in ([5]) we have o^Uo = Ao- Using the non-singularity of the matrices 
o.q and Ao, IXo is a non-singular matrix. In an analog way, taking m = 1 in (J5|), we 
have 

ckqIXo + a\ IXi = Onxn ■ n, lr -iir ^ A 

alu 1+ alU 2 = A u ^ - U X U U x ) = A x . 

Since Ai and a\ are non-singular matrices then det(!X2 — IXiIXq^IXi) 7^ and, as a 
consequence, the second leading principal submatrix is non-singular. This argument 
can be inductively used and we obtain that IX is quasi-definite. 
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Conversely, to find the vector sequence of polynomials such that {'B m } m6 N with 
23m = Sj=o a T'^^ where aj 1 <E Mat x at(R) and where a™ is non-singular lower 
triangular matrix such that (h k U) ( r B rn ) = A m 4, m , k — 0,1,..., m, m € N, is 
equivalent to solve 

" u(iPo) ■•• u(y m ) 

[ < a? ■■■ a™ ] : ••. j 

_ U(7 m ) • • • U(7 2m ) 



[ 



For m = 0, we have a^Uo = Aq. Using the non-singularity of Uq, and the decom- 
position LU , we can find uniquely a® a non-singular lower triangular matrix and 
Ao a non-singular upper triangular matrix such that o^Uo = Ao. 

For m = 1 we have 

i i„ , a i.e. ai(u 2 -UiU Ui) = Ai. 

ag Hi + U.2 = Ai, 1V u 

Again, using that the second leading principal submatrix U2 — lliUg 1 lXi is non- 
singular and the LU decomposition we can find uniquely a\ a non-singular lower 
triangular matrix and Ai a non-singular upper triangular matrix such that a\ — 
(U2 — UiUq 1 Ui) — Ai. We also obtain from a\ Ho + a\ Ui = 0at x tv, uniquely 
the matrix aj. This argument can be inductively used and we obtain the stated 
result. □ 

Theorem 4. Let U be a quasi- definite vector of linear junctionals and let {23 m } mS N 
be a sequence of vector polynomials. Then, the following statements are equivalent: 

(a) The vector sequence of polynomials {23 m } me N is left- orthogonal with respect 
to the vector of linear functionals IX, i.e., 

(h k U) (B m ) = A m S k , m , k = 0,1,. . . ,m, m G N, (10) 

with A m a non-singular N x N upper triangular matrix given by 

A m = C m ■ ■ ■ Ci A , m > 1, 

where Ao is a N x N non-singular matrix and {C m }meN is a sequence of 
non-singular upper triangular matrices. 

(b) There exist sequences of N x N matrices, by {A m } m6 pj,{.B m } me N, an< ^ 
{C m }meN with C m a non-singular upper triangular matrix, such that 

h(x)-B m (x) = A m %n+l{x) + BmSmix) + C m %n-l{x), 771 > 1, (11) 

with = 0at x i and 23o(x) = 7q(x), where 7q(x) — [l x ■ ■ ■ i^ 1 ] . 

Proof. To prove (a) =>■ (6), first we consider the polynomial hT> m and then we take 
into account that hT> m is a polynomial of degree m + 1 that can be written 

m+l 

Ma;)3 TO (aO - 51 A k^k(x), A™ e Mat x jv(R). (12) 

fc=0 

We will prove that A™ = 0at x jv, for fc = 0, 1, . . . , m — 2. Indeed, if we apply the 
vector of functionals U to both sides of (fl~2|) then we get A™ — 0n x n- 
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Thus we can rewrite (fT2|) 

m+l 



h(x)'B m {x) = Y, A k^k{x). (13) 

k=l 

Again, by applying the vector of functional U to both sides of (fT3]) we obtain 

A™ = Onxn- 

Iterating this procedure, i.e., first by multiplying by h 2 , afterwards by h 3 and then, 
successively, applying U we get 

A™ = NxN , for fc = 0, l,...,m-2. 

Thus, we can rewrite (|12l) 

h(x)-B m (x) = AZ^m-xix) + A™¥> m (x) + A" +1 £ m+1 (x). (14) 

If we multiply (fl4|) by /i" 1-1 and we apply the vector of linear functionals U then 
we obtain 

i4£-i = (^) (Sm) ((/i^'lX) (S,^!))" 1 = A ro A- 1 _ 1 m > 1 . 
Using the same technique we obtain 

AZ = [(/i m+1 l[)(a3 m )-A™_ 1 (^»U)(15 m _ 1 )] [(^U)(S m )]- x 
= [(^ +1 1X) (2 m ) - A m A m 1 _ 1 (ft™U) (S m _ x )] A- 1 

and - [(/i™+ 2 U) (B m ) - A™_ 1 (/ l m+1 U) - A™(/ 1 " l+1 lt) (S m )] A^. 

The comparison with the coefficients in (jlip yields the following explicit expressions 
for the coefficients in the recurrence relation: 

A m = [(h m+2 U) (S m ) - AZ^(h m+1 U) (B m _x) - A™(/, m+1 U) (S m )] A m x +1 , 
B m = [(ft m+1, U) (25 m ) — A™_ 1 (/i m U) (B m _i)] A^ 1 , 
Cm = A m A m _-^. 

To prove (b) (a), we must start by constructing a vector of linear functionals U, 
satisfying (TTOj) . which is defined from the sequence of moments {U m }meN using 

U(S ) = Ao, U(B TO ) = M , m> 1, (15) 
where Ao is a non-singular upper triangular matrix. Since {"Pj}jeN with 

^-(.t) = (h(x)yy a {x) and 3» (a;) = [is • ■ • i" _1 ] T , 

is a basis for P^, then there exists a unique family of matrices 7™ € Mtv x at (R) 
such that the vector of polynomials CB m can be written "B m (x) — Y^jLo lj l '^'j( x )- 
Thus, 

. For m = 0, U(B ) - To^o), i.e., U = (rg^Ao. 
1 

. For m = 1, U(Si) = X^-U (?,-), i.e., U x = -(Ti^To^o- 
3=0 

2 1 
•Form = 2, U(S a ) = ^Ti^i), i.e., U 2 = - ^(tI) -1 ^- 
3=0 j=o 
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m— 1 

For m > 3, we have ll m = - ^ (7™)~Sf ^ . 

j=o 

First, we will prove that for II defined as above, we have 
{h k U) (B m ) = NxN , m>k+l. 
To prove it, we apply II in the recurrence relation: 

U [hfB m ) = A m U (S m+1 ) + B m U (S m ) + C m ll (B m _i) = M , m > 2 . 
Again, if we multiply both sides of the recurrence relation by h, then we obtain 

h 2 (x)?i m (x) = h(x)A m 'B m+1 (x) + h(x)B m 'B m (x) + h(x)C m 3 m -i(x), 
and, as a consequence, by applying IX in the last relation, we have for m > 3 

(h 2 U) (S m ) = A m (MI) (S m+ i) + B m (/ill) (B m ) + C m (MI) (S m _i) = O^vxw . 
Proceeding in a similar way we get 

(h k U) (S TO ) = Oat x7V , m > fc + 1, i.e., (/i fc U) (S TO ) = O^xat , fc = 0, 1, . . . , m - 1. 
For k = m, we have 

(h m U) (% n ) = A^h^U) (S m+1 ) + B^h^U) (S m ) + C m {h m - l U) (B m _i) , 
and so 

(^ m U) (S ro ) = C m (ft ro - 1 U) (S m _ x ) = C m C TO _ x • • • CiAo , m > 1 . 

Therefore, the moments associated with the vector of linear functionals IX are 
uniquely determined from (|15p . Thus, we obtain the orthogonality conditions (|10p . 
Hence, the result follows. □ 

Next we will introduce the concept of right-orthogonality with respect to a vector 
of linear functionals and, afterwards, we will show how the right and left vector 
orthogonality are connected. 

Definition 3. Let II = [it 1 • • • u N ~\ T be a vector of linear functionals and let 
consider a sequence of matrix polynomials {G m } m £N- {G m } m gN is said to be right- 
orthogonal with respect to the vector of linear functionals IX if 

(a) deg G m = m. 

(b) (Gl(h(x))U) (Vj) = NxN , j = 0, 1, . . . ,m - 1. 

(c) (Gf n (h(x))1i) (T TO ) = 6 m , m e N, where 9 m is a non-singular lower trian- 
gular matrix. 

Concerning the right-orthogonality we obtain some analog results to those we 
found for left-orthogonality. For example, the matrix right-orthogonal polynomial 
sequence is uniquely defined up to a multiplication on the right by a non-singular 
matrix and the vector of linear functionals IX is quasi-definite with respect to the 
right-orthogonality. We will present these results but we shall skip the proofs since 
the techniques are the same as used in the left-orthogonality case. 

Theorem 5. Let 11 be a vector of linear functionals. Then, IX is quasi- definite if 
and only if there exists a sequence matrix polynomials {G m } mS N, with G m {h{x)) = 
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Sjlo ftj 1 '(h-( x )y ') f or Pj 1, G 3VCatxA'(R) where (3™ is non-singular upper triangu- 
lar matrix and there exists a sequence of non-singular lower triangular matrices, 
{6m}meNi such that 



Moreover, 



G m =[l h 





= 0, 1, . . . , m - 


"I, 


m G N. 


U(To) • 




-1 





U(? ro ) • 


U(J l 2 m ) 




. ™ . 



(16) 



Theorem 6. Let 11 be a vector of linear junctionals and let {G m } mC N be a sequence 
of matrix polynomials. Then, the following statements are equivalent: 

(a) The sequence of matrix polynomials {G m } m eN is right- orthogonal with re- 
spect to the vector of linear junctionals U, i.e., 

(GZ(h(x))U) (Tj) = NxN , j = 0,l,...,m-l, 



(Gl(h(x))U) (0> m ) 



e, 



to G N, 



where G m is a non-singular lower triangular matrix. 
(b) There exist sequences of N x N matrices, {D m } me N, {-Em}meN> and 
{-Pm}meN wi/i -Fm a non-singular lower triangular matrix, such that 

h{x)G m {h(x)) = G m+ i(h(x))D m + G m (h(x))E m + G m -i(h(x))F m , m > 1 , 

im'i/i = On xN and Gq(x) = Inxn ■ 

To show the connection between right and left vector orthogonality we will 
introduce some concepts on duality theory. We denote by P* the dual space of P, 
i.e., the vector space of complex valued linear functionals defined on P. 

Let {p m }meN be a sequence of scalar monic polynomials. The sequence of linear 
functionals {i/ n }neN, where L n G P* is said to be its dual sequence if L n (p m ) — 
S m ,m Tn, n G N, where 5 n ^ m is the Kronecker delta. 

Let {L n } n( £fi be a sequence of linear functionals. The vector sequence of linear 
functionals {£ ra }neN given by 

£n = [L n N ■ ■ ■ i( n +l)JV-l] i n G N , 
is said to be the vector sequence of linear functionals associated with {L n } n( z^. 
Taking into account Definition [TJ we get 

LnNKPmN) "" Lr n+1 - )N _ 1 (p mN ) 



NxNOm.n ■ 



L n N{P{m+l)N -l) ' ' ' £(n+l)A'-l(??(m+l)A'-l)_ 

Definition 4. Let {23 m } me N be a vector sequence of polynomials. The vector 
sequence of linear functionals {<C n }neN i s sa id to be its dual vector sequence if 

UnCBrn) = iNxN $m,n, U, TO G N . 

The next two results give the connection between right and left vector orthogo- 
nality through the equivalent conditions of these two types of vector orthogonality. 
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Definition 5. Let IX be a vector of linear functional. We denote by It, the nor- 
malized vector of linear Junctionals associated with IX, IX = ((1I(!Po)) _1 ) T 'U , where 
? (x) = [l,x, ■■■ ,x JV - 1 ] T . 

Furthermore, from this definition we have 

U(3> ) = (((U(J'o))' 1 ) T U)(J'o) = UCToJWTo))" 1 = Inxn- 

Theorem 7. Let IX be a quasi- definite vector of linear functionals, {!B m } me s be a 
sequence of vector polynomials, and let {fj n }n&i be its dual vector sequence. Then, 
the following statements are equivalent: 

(a) {S m } me N is left- orthogonal with respect to IX. 

(b) There exist sequences of N x N matrices, {^4 m } m eN, {-Bm}meN, and 
{C m } me N with C m a non-singular upper triangular matrix, such that 
{23m}meN satisfies the three-term recurrence relation 

h{x)'B m {x) = A m T, m+1 {x) + B m <E, m {x) + C m S m _i(x), m > 1 (17) 

with = Oi x tv and Ho(x) — fo(x) , where fo(x) — [l x ■ ■ ■ x N ^ x ] . 

(c) There exist sequences of N x N matrices, {A n }^, {-B n } n6 N, and {C n }„ g N, 
with C„+i a non-singular matrix, such that {£j n }n£N is defined by the three- 
term recurrence relation 

h{x)L n = {C n+1 ) T L n+1 + {B n ) T L n + {A n _ l ) T L n _ l , n > 1, (18) 
where L = ((U(9 )) T )- 1 U, £ x = (Cf)- 1 ^)' - (B ) T )[IX(5' )) T ]- 1 U . 

n 

(d) There exist matrix polynomials G n (h(x)) with G n (h(x)) — / ^ (3^ (h(x)y , 

3=0 

where /3™ is a non-singular matrix, such that the elements in the dual vector 
basis {XL n } ra gN can be written in terms of the vector of linear functionals IX 
as follows 

Hn = {G n (h(x))f U, neN. (19) 

(c) The sequence of matrix polynomials {G n }„gN defined by (|19|) satisfies 

h(x)G n (h(x)) = G n -i(h(x))An-i + G n (h(x))B n + G n+1 (h(x))C n+1 , n > 1 , (20) 

loz'tft, initial conditions G-i(h(x)) — Onxn and Go(h(x)) — IX(I'o) 1 ■ 
(f) The sequence of matrix polynomials {GnJneN defined by (|19|) is right- 
orthogonal with respect to the normalized vector of linear functionals IX. 

Proof. We will prove this theorem according to the following scheme: 

(a) (6), (e) «■ (/), (6) ^ (c), (c) (d), (d) => (e), and (e) => (c). 

The proofs of (a) O- (6) and (e) 4=> (/) follow immediately from Theorems 0] 
and [5] We start by proving that (b) (c). Let 

n+l 

h(x)L n = Pj^i, where (fff = = M^)%) = J € N. 

j'=o 
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Applying the vector of linear functionals £„ in both sides of the three-term recur- 
rence relation satisfied by {Hk}keN, we have 

(/3™) T = A k L n CB k+1 ) + B k L n (V k ) + C k L n (V k ^)) 

An-i, j = n-l, 

B n , j = n, 

C n +i, j = n+l, 

Ojvxjv, j ^ n- l,n,n+ 1, 

i.e., = A^_ 1: = Bj x , and = Cj +1 . Thus the three- term recurrence 

relation for the vector sequence of linear functionals {-C„}ngN follows. 

m+l 

To prove that (c) => (b), let h"B m = ^ 7™% 7™ 6 M NxN (R). Applying 

j=o 

the vector linear functional L n in both sides of the last relation, we get /iiL„(23 m ) = 

rn+l 

7™,C n ('Bj) = 7™. Now, from our hypotheses, we have 

lrt = f*>n+\\&m)C'n+\ + h n ('& r n)B n + £ n _i)(!B m )j4„_i 

G m , ra = m - 1, 

5 n , n = m, 

A m , n = m + l, 

Onxn, n^m-l,m,m + l. 

So, ([IT]) holds. 

To prove (c) =>- (<f) we will show by induction that {£n}neN has the follow- 
ing representation L n = (G n (h(x))) 11, n G N. For n = 0, we have that Lq = 
((U(CPo)) T ) _1 lI. Now, let us assume that the statement holds for k = 0, 1, .. . ,p, 
i.e., Z k = (G k (h(x))) U with degG k = k, k = 1,. .. ,p. We will show that it is also 
true for k — p+1, i.e., L p+ i = (G p+ i(h(x))) T U, p £ N. Considering the three-term 
recurrence relation satisfied by {iL p } p gN and taking into account the hypothesis of 
induction, we have 

£ P+1 = (C P+ i)- T [(h(x)I - (BT))G p (h(x)) - (A p ^) T G p ^(h(x))] U 
= [(G p (M^)) T ((^(^K-(i3 P ))-Gj_ 1 (/ i (x))^ 1 ))C p+ 1 1 ] T U. 

Thus, = {G p j r i(h(x))) T It, p e N, i.e., if the condition holds for fc = 1,. .. ,p, 

then it is also true for p + 1 . 

To prove that (d) =4> (e) we will write hG^ in terms of {Gj} je N, i-e., 

n+l 

h{x)Gl{h{x)) = Y,^G^{h{x)), where a™ e MAr xA r(M). (21) 
Thus, the multiplication on the right by II in both sides of the last equation yields 

n+l 

h(x)Gl{h{x))U = ^2a^Gj(h(x))U. 
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Applying this relation to 23 k we get 

n+l 

{h{x)G T n {h{x))U)(<2> k ) = £(a7GT(M»))U)(S fc )- 



3=0 

Since, £.„ = (G„(/i(:r))) T U, neN, 

n+l 

L n (/i(z)23 fe ) = J2 ^ K) T = «) T 
i=o 

Using (JTTJ) in (JUl) we get 

K) T = C k L n {'B k - 1 )+B k L n {'B k )+A k Z n {'B k+1 ) 

A n —i-, k = n — 1, 

B n , k ~ n, 

Cn+i, fe = n + l, 

OwxAf, fc^n— l,n, n+l. 

Thus, {G„} nS N satisfies 

ft,(a;)G„(/i(ar)) = G„_ x (/i(x))A„_i + G n {h{x))B n + G n+1 (h(x))C n+1 . 

Finally, to prove that (e) => (c), we must take the transpose in the recurrence 
relation (f20)) and, then, multiply on the right by II both sides of the resulting 
equation. Thus (|18p follows. □ 



Theorem 8. Let 11 be a quasi- definite vector of linear junctionals, {23 m } me N be a 
sequence of vector "polynomials and {G n }„ e N defined by flB]) . Then, { l B m }meN and 
{G„}„ e N are bi-orthogonal with respect to 11, i.e., 

{{G n {h{x))) T U){% n ) = I NxN S n>m ,n, to e N, 

if and only if A m = otic? f2„ = (a™) -1 . 

As a consequence, the dual sequence {£i n }„gn associated with {25 m } mg N is given by 
&n = {G n (h(x))) T U, neN. 

Proof. There exists a unique family of matrices (a™) C Mat x at(R) such that 23 m — 
Y^ijLo a j l ^'ji where a™ is a non-singular matrix. Hence, 

m m 

(G^(M*))U)(S m ) = {G T n {h{x))U)(£af'S> j ) = £af (G^))^?;)- 

3=0 3=0 

Since {G n } ne N is right-orthogonal with respect to IX then 

{G T n {h{x))U){'B m ) = { < 0m ' "* = " 

Thus, {G^ n {h{x))U){ r B m ) = I NxN if and only if a™ Q m = InxN, i.e., 9 m = (a™)" 1 . 
Now, let us consider 

n n 

(GZ(h(x))U)(B m ) = (C^PjiHxWfU)^) = ^(WlX)(S ro )^. 

3=0 j=0 
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As above, since {'B m } me N is left-orthogonal with respect to IX then 

(GZ(h(x))U)(V m ) = \ ™ = n 

L UjvxAf, m > n. 

So, (Gl(h(x))U)(B m ) = I NxN if and only if A m = □ 

To conclude this section notice that Theorems [2] and [7] suggest that the se- 
quence of matrix polynomials {y m } m eN is orthogonal with respect to some matrix 
of measures. As a consequence of Theorem [3 the sequence of matrix polynomials 
{G m } mS N should also be orthogonal in the matrix sense. Finally, Theorem [8] sug- 
gests that the sequences of matrix polynomials {V m } m( z^ and {G m } m(E N should be 
bi-orthogonal to each order, as we will prove in the next section. 

3. The connection between vector and matrix orthogonality 

In this section we show how the vector and matrix orthogonality are connected 
when a special case of the matrix orthogonality is considered. In this sense, if 
the sequences {G n }„ e N and {23 m } mS N are bi-orthogonal with respect to the vector 
of linear functionals IX, then the sequences of matrix polynomials {G ra } ng N and 
{KnjmeN are bi-orthogonal with respect to a complex matrix of measures. 

Definition 6. Let It be a vector of linear functionals. We define the generalized 
Markov matrix function, 3", associated with IX by 



J(z) := U x 



y (x) 

z — h(x) 



1 ^ ■-- (u? 1 



z-h(x)l \ x ) z -h{x)i 



(22) 



N***) z-h(x)l \ x ) z -h{x)l 

with z such that \h(x)\ < \z\ for every idL where L = U^i,...^ suppttj . Here U x 



represents the action of IX on the variable x and To (a;) = [1 x 



x N - 1 } 



Theorem 9. Let U be a quasi- definite vector of linear functionals and let 5" be its 
generalized Markov matrix function. Then, the following statements are equivalent: 

a) The sequences {G n } n gN and {B m } me g are bi-orthogonal with respect to IX, i.e. 

((G n (h(x))) T U x )(B m ) = I NxN 5 n>m , n,m e N. 

b) The sequences {G n } ne n and {V m } me ®, where ¥> m (z) = V m {h{z))'?o(z) , are 
bi-orthogonal with respect to J, i.e., 

7T— / V m (z)'J(z)G n (z)dz = I NxN 6 n ,m, n, m € N. 
2m J c 

where C is a closed path in {z G C : \z\ > \h(x)\,x S L}. 
Proof. Taking into account that 
V m{ z)H Z )G n (z) = (^^) U-G n ( Z ) = ((G n (z)fU x ) (^g^ 

l - jj m {z)S{z)G n {z)dz = ^-.j c {{G n {z)) T U x ) ( ^yff 



we have 
1 

~2m 



dz . 
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Because of G n , V m , and To are analytic functions, according to the Cauchy integral 
formula we have 

-L / ((G n (z)) T U x ) ( V ^ z )^) \ dz = ((G n (h(x))) T U x )(V m (Mx))y (x)), 
atti j c \ z — n(x) j 

and, as a consequence, for all n, m € N 

-i- / V m (z)9(z)G n (z)dz = ((G n (h(x))fU x )CB m (x)) = I NxN S n , m . 
Thus the statement follows. □ 



The last theorem tell us that {CB m }„ l6 N is a sequence of vector polynomials left- 
orthogonal with respect to U if and only if {V m } m £fi associated with {'B m }rneN is 
left- orthogonal with respect to 3\ Also, {G m } me N is a sequence matrix polynomials 
right- orthogonal with respect to U if and only if {G m } me N is right- orthogonal with 
respect to SF. 

It is important to recall now that the definition of J shows us that we only need N 
linear functionals to describe the matrix orthogonality. Usually, to describe the 
matrix orthogonality, (1 + N)N/2 measures are needed (see, for example [3JH1IS]). 

As we have already referred in the introduction, we need to know when a se- 
quence of matrix polynomials defined by a recurrence relation ([7]) is related to some 
kind of matrix orthogonality. Partial answers were given to this problem, but no 
complete answer was given as far as we know. To do that, we start by considering 
a N x N matrix of measures W that is not necessarily positive definite in R, and 
such that there exist matrix sequences {V m } m( zfi and {G m } mg N, orthogonal with 
respect to W in the following sense 

/ V m (x)dW(x)x k = n^5 k , m , fc,m>0, (23) 
J x k dW(x)G m (x) = r» 2 m 4,m, fc,m>0, (24) 

where f^ n is a non-singular upper triangular matrix, Q 2 m is a non-singular lower 
triangular matrix, and <5/c,m is the Kronecker delta. 

V m and G m are matrix polynomials of degree m with non-singular leading co- 
efficients and they are defined up to the multiplication on the left or on the right 
by a unitary matrix, respectively. The matrix sequences {Vm} m eN (respectively, 
{G m } m6 N) satisfying (|2"3"]) (respectively, (25)) are said to be the left- orthogonal 
matrix polynomial sequence (respectively, right-orthogonal matrix polynomial se- 
quence), with respect to the matrix of measures W. 

Usually, in the theory of matrix orthogonal polynomials there are only references 
to the left-orthogonality. The reason is that the authors deal only with orthonor- 
mality with respect to a positive definite matrix of measures, that allows us to say 
that left and right orthogonality are, essentially, the same. Very few authors have 
emphasized this difference (see, for instance, [TU]V 

The moments of the matrix measure W are given by N x TV matrices 
S k = / x k dW{x), k = 0,l,.... 
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From the orthogonality conditions it follows that the sequences {V m }meN and 
{G m }meN satisfy three-term matrix recurrence relations. It is not so obvious to 
prove the converse result, i.e. if a sequence of matrix polynomials is defined by a 
recurrence relation (J7J or (|20|) . then there exists a matrix of measures W, not nec- 
essarily positive definite, such this sequence is left or right-orthogonal with respect 
to W. 

The following result proves this equivalence with respect to the left-orthogonality 
and gives an extension of the Favard's theorem in the matrix case. 

Theorem 10. Let {V m } m <aN be a sequence of matrix polynomials. Then, the fol- 
lowing statements are equivalent: 

(a) The sequence {Vm} m eN is left- orthogonal with respect to matrix of mea- 
sures W . 

(b) There are sequences of scalar matrices {A m } me ^, {B m } me ^ and {C m } me ^, 
with A rn lower-triangular, and C m+ i upper-triangular, non-singular matri- 
ces for 
m G N, such that the sequence {V m } m ^fi satisfies 

zV m (z) =A m V m+1 {z) +B m V m (z) + C m V m - 1 (z), m > 1, (25) 

where V-i(z) = OnxN and V$(z) = InxN ■ 

Proof. First we will prove that (a) implies (b). Since the sequence {Kn} m6 N is a 
basis in the linear space of matrix polynomials we can write 

m+l 

zV m (z) = A™V k (z), A™ e M NxN (R). 

fc=0 

Then, from the orthogonality conditions, we get 

A T [ V 3 {z)dW{z)z j = [ V m (z)dW(z)z j+1 = NxN for j = 0, . . . ,m - 2. 
Thus, 



zV m (z) = AZ^V m -i{z) + AZV m (z) + A™ +1 V m+1 (z), 

^zV m dW(z)z m ^j (fi^)" 1 , A™^ = (^J zV m (z)dW(z)z m - 

x(nl^)-\ and AZ +l = ^zV m (z)dW(z)z m+1 ^ (Qi, +1 ) _1 . 

Taking A m = A% +1 , B m = A™, and C m = A™_ x the result follows. 

Finally, to prove that (6) implies (a), we should start by defining recursively 
the matrix moments associated with the matrix of measures W by the following 
conditions 



Sq 



[ dW{z) = nl and / V m (z)dW(z) = Onxn, m>l, 



where fij is a non-singular upper triangular matrix. 
Taking into account that V m can be written 

V m (z) = V m ,mZ m H h V m .iZ + Vmfi 
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with V m ,m a non-singular matrix, then we have 

OnxN = / V m (z)dW(z) = V m>m Sm H 1- VmflSo- 



Thus, the moments are defined recursively by S m — V m ^ m V m jSj . 

3=0 

Let us show that 

V m (z)dW(z)z k = Oatxjv, fc = 0,...,m-l, and / V m (z)dW(z)z m = Q^. 



From we get / zV m (z)dW(z) — Onxn, to > 2, Again, by multiplying both 
sides of the recurrence relation by z we get 

z 2 V m (z) = A m zV m+1 (z) + B m zV m (z) + C m zV m ^ 1 (z), 
and , as a consequence, J V m (z)dW(z)z' 2 — Onxn, to > 3 . 

In an analog way, we conclude that / V m (z)dW(z)z k — Omxn, to > fc + 1, and 

Jr. 

so / V m (z)dW(z)z k = NxN , k = 0,...,m-l. 

JR 

For fc = m we have 

V m (z)dW(z)z m = C m [ V m ^(z)dW{z)z m - 1 =C m C m . 1 ---C 1 nl 

I JR 

Then 

V m {z)dW(z)z k = NxN , fc = 0, . . . , to - 1 and f V m (z)dW(z)z m = 9} m , 

: JR 

where 57^ = C m C TO _i • ■ ■ CifiJ is a non-singular upper triangular matrix. □ 

As an important remark the reader should notice that left vector orthogonality 
and matrix orthogonality are equivalent. This equivalence is given by Theorems[21[7l 
and 1101 A similar result can be obtained for the right-orthogonality. 

Theorem 11. Let {G m } m eN be a sequence of matrix polynomials. Then, the fol- 
lowing statements are equivalent: 

(a) {G m } me N is a right- orthogonal sequence of matrix polynomials with respect 
to a matrix of measures W . 

(b) There are sequences of scalar matrices {A m } m( z^, {-B m } m eN and {C m } m£ N, 
with A m -i lower-triangular, and C m +i upper-triangular, non- singular ma- 
trices for 
m £ M, such that the sequence {G m } m gN satisfies 

G m {z) = G m -i(z)A m -i + G m {z)B m + G m+1 (z)C m+ i, m>l (26) 
where G-±(z) = Onxn and Gq(z) — Inxn ■ 
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4. Some characterizations of the vector and matrix orthogonality 

In this section we present some characterizations of the vector orthogonality 
as well as of the matrix orthogonality. First, we analyze two type Hermite-Pade 
approximation problems that characterize completely the right and left vector or- 
thogonality, respectively. 

Definition 7. Let {!B m } m£ N be a vector sequence of polynomials, let {G m }mefi 
be a sequence of matrix polynomials, and let U be a quasi-definite vector of linear 
functionals. The sequence of polynomials {23^} me pj given by 

:=U. (M^f„(.) 

is said to be the sequence of associated polynomials of the first kind for {'B m }meN 
and U. In a similar way, the sequence of polynomials {Gm}meN given by 



Grn+l( z ) _ G m+l {h(x)) 



CM*)), 



z — h(x) 

is said to be the sequence of associated polynomials of the first kind for {G m } m gN 
and U. Here U x represents the action of U on the variable x. 

Theorem 12. Let XL be a quasi- definite vector of linear functionals, {1>rn\m&H a 
vector sequence of polynomials, {Sm } m£ j its sequence of associated polynomials of 
the first kind, and 3 the generalized Markov function given in (|22p . Then {'B m } mg f!i 
is left- orthogonal with respect to the vector of linear functionals li if and only if 



V m+ i(z)J(z) ~ 03«(z) = A m+1 -^ + 



Proof. From the definition of B^' , we get 

B &>(*) = U x ( = Vm+l{mz) _ Ux 

V z - K x ) J 

But 

"B m+1 {x)\_^ (^{h{x)T m f _ s \ _ ^ ((h( x)rU x )(B m+1 (x)) 

r n+l 



z — /i(x) 



Hence, 



/^+i(x)\_ ~ ((fe( a: ))"lX a .)(23 m+1 (x)) _ 1 



7 n— JTfc+1 

if and only if the sequence {23 m }meN is left-orthogonal with respect to 1i x . □ 

Theorem 13. Let li be a quasi- definite vector of linear functionals, {G m } mg N 
a sequence of matrix polynomials with N x N matrix coefficients, {Gm } m£ n its 
sequence of associated polynomials of the first kind, and SF is the generalized Markov 
function. Then {G m } m6 N is right- orthogonal with respect to the vector of linear 
functionals li if and only if 



J(z)G m+1 (z) - G«(z) = 6 m+1: ^ + 
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Proof. Taking into account the definition of the polynomial {Gm } me s, we have 



G&Hz) = 



GL+i(z)-C^ +1 (h(x)) 



U x 



z — h{x) 

u * {-^Wj) Gm+l[z) (- W 



But, 



Hence, 



z — h(x) 



(G T m+1 (h(x))U x ) (^y) =e ro+1 ^ + ••• , 
if and only if the sequence {G m } mfE N is right-orthogonal with respect to VL X . □ 

Next, some algebraic results concerning the behavior of the sequences of matrix 
orthogonal polynomials {V m } me ^ and {G m } me N are given. 

Theorem 14. Let h be a polynomial of fixed degree N and U be a quasi- definite vec- 
tor of linear functional^. Let {G m } me N d/i<i {Bm}men be, respectively, sequences of 
matrix polynomials with deg G m — m, for all m € N and ¥> m (x) — V m (h(x))'J 3 o(x) , 
where V m is a matrix polynomial with degV^ = m, for all m e N. Then, the 
following statements are equivalent: 

(a) {23 m }„ l£ N is a sequence of vector polynomials left-orthogonal with respect 
to It. 

(b) {X„} rie N is a sequence of vector linear functionals bi-orthogonal with respect 
to {23 m } m( =N such that L n = G^(h(x))11. 

(c) {V^j} m£ N and {G m } me N satisfy the Christoffel-Darboux type formula 

m 

{x - z)^G k {z)V k (x) = G m (z)A m V m+1 (x) - G m+1 (z)C m+1 V m (x), (27) 

fc=0 

with x, z e C. 

(d) For every m e N, {Kn}meN an d {G m } m£ N satisfy the confluent formula 

~ G m+ i (x)G m+ iy m (a;) = W xAf, (28) 

m 

^Gfe^t/fe^) = G m (x)A m y: +1 (x) - G m+ i(z)C m +i^m(*), (29) 

fe=0 

(e) For every m e N, {Kn}meN a^d {G m } m£ N satisfy for all ifC 

G m (.T)A m y m+ i(a;) - G m+ i(a;)G m+ iV r m (a;) = NxN , (30) 

£Gt(i)V t (i) = G' m+1 (x)C rn+1 V m (x) - G' m (x)A m V m+ i(x) . (31) 

k=0 

Proof. To prove this theorem we will proceed according to the following scheme 
(a) & (6), (b) => (c) (e), (e) => (b) and (c) (d) => (a). 
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The equivalence (a) (b) is proved in Theorem [71 To prove that (b) implies (c) 
we remember that the sequences of matrix polynomials {V m }meN and {G m }men 
verify, respectively, the recurrence relations 

xV m {x) = A m V m+ i(x) + B m V m (x) + C m V m -i(x) (32) 
zG m (z) = G m -i(z)A m - 1 + G m {z)B m + G m+ i(z)C m +i (33) 

Multiplying on the left by G m (z) in both sides of (f3"2"]) and on the right by V m (x) 
in both sides of (j3"3")l and subtracting the resulting expressions, we get 

(x - z)G m (z)V m (x) = [G m (z)A m V m+ i(x) - G m - 1 (z)A m ^ 1 V m {x)] 

- [G m+ i{z)C m+1 V m {x) - G m {z)C m V m -x(x)} 

and so we have (|27P- To prove that (c) implies (d), we just have to take z = x in (|27|) 
and then we obtain (f2"5| . (|29|1 follows from (|2"T|) by differentiation with respect to 
x and letting z = x. 

To prove that (c) implies (e), we must take z — x in (|2T[) and then ((30)) holds. (|3Tj) 
follows in a similar way by differentiating (|2T[) with respect to z and taking z = x. 

To complete the proof we need to show that (d) implies (a) . We can rewrite the 
equation (j2T))) as 

m— 1 

G m (x)A m V^ +1 {x) - G m+1 {x)C m+1 V^(x) = G m (x)V m (x) + ^ G k (x)V k (x), 

or, equivalently, 

G rn (x)V m (x) = G m (x)[A m ^ l+1 (x) + CmV^ix)} 

- [G m +i(x)C m +i + G m -i(x)A m -i]Vj n (x). 

Using (HH1) we get [(A m F m+ i(a;) + C m V m -i(x))V m 1 (x)}' = I N xN ■ Then, we have 

[A m V m+1 {x) + CmVm-^x^V- 1 ^) =xl- B n , 

i.e, {V m } m £fi satisfies a three-term recurrence relation. Now, multiplying both 
sides in the three-term recurrence relation by CPo, from the definition of 23 m and by 
Theorem SI the result follows. 

Finally, to prove that (e) (b) we proceed in a similar way as in the proof of 
(d) (a) starting from (|3Tj) and taking into account Theorem [Jj □ 

5. Markov type theorem 

The block matrix 

Bq A Onxn 
Ci B l A l 
Onxn C2 B 2 

is related to the matrix polynomial sequences {V m } me ?i and {G m } m£ N trough the 
recurrence relations (J7J and (j2^]) . This block matrix is said to be the iV-block 
Jacobi matrix associated with the above matrix polynomial sequences. 
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For polynomials satisfying a symmetric recurrence relation, it was proved in [8] 
that the zeros of the TO-th matrix orthogonal polynomial are the eigenvalues of the 
leading principal submatrix J m N of J. This result can be generalized for sequences 
of orthogonal polynomials that satisfy non-symmetric recurrence relations. Thus, 
for m S N, the zeros of the matrix polynomials G m and V m are the eigenvalues 
of the matrix J m N (with the same order of algebraic multiplicity) where I m NxmN 
is the mN x mN identity matrix and J m N is the leading principal submatrix of 
dimension mN x mN for the iV-block Jacobi matrix. 

Lemma 1. 5J Let V{t) be a N x N matrix polynomial and let a be a zero ofV(t) 
of multiplicity p, i.e, a zero of multiplicity p of the scalar polynomial detV(t). Let 

L{a,V) = {v € C N : vV{a) = Q lxN } and R(a, A) = {v e C N : V(a)v* = N xi}- 

If dim L(a, V) = dim R(a, V) — p, then 

(Adj (V(t))f } (a) = NxN , for I = {),... ,p-2 and, (Adj (F(t))) (p ~ 1} (a) ^ NxN . 

Moreover, rank (Adj {V(t))) l ' p ~ 1 ' > (a) = p and (Adj (l / (t)))' p-1 - ) (a) defines a linear 
mapping from C N onto L(a, V) which is an isomorphism from R(a, V) into L(a, V). 

Lemma 2. Let x m _k, k — 1, . . . , s with s < mN be the zeros of the matrix poly- 
nomial V m . For any matrix polynomial V(t) of degree < n — 1 we have the partial 
fraction decomposition t £ C \ {x m ,i, ■ ■ ■ , %m,s}i 

K—l 

where C m k = ^ -V"(x m . fc )(Adj (V m (t)))( lk ~V (x m k ) and l k is 

(det(F m (t)))t'*J(ai m ,fc) 

the multiplicity of x m> k, (lk < N). 

With these results we are able to establish a quadrature formula for the matrix 
orthogonal sequence {KnjmeN- 

Theorem 15 (Quadrature Formula). Let {Kn}meN oe the sequence of matrix poly- 
nomials that is left- orthogonal with respect to the matrix of measures W . Also let 
{'S> m }men be the sequence of vector polynomials defined by ((4]) and let {Dm } ra eN, 
be the sequence of associated polynomials of the first kind for {'B m } me N and IX. Let 
%m,k, (k = 1, . • • , s) be the zeros of the matrix polynomial V m (hence s < mN), and 
let r m> fc be the matrices 

Tm - k = I a Tn? < ^ {V m {x)))^- 1 \x m . k )Ti^l 1 {x m , k ) , 

(det (V m (x))y l *>{x m . k ) 

for k — 1, . . . , s where l k is the multiplicity x m . k . 

Then, for any polynomial V of degree less than or equal to 2m — 1 the following 
quadrature formula holds 

s 

V(h(x))dW(h{x)) = J2 V{x m , k )T m>k . 

k=l 
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Proof. Let V be a matrix polynomial of degree less than or equal to 2m — 1 . Since 
V m is a polynomial with non-singular leading coefficient, then (cf. llj) 

V(x)=C(x)V m (x)+R{x), 

where C and R are matrix polynomials with degree of R less than or equal to m—1. 
Thus V(x)V~ 1 (x) — C (x) + R(x)V r ^ 1 (x) assuming that x is not a zero of V m . Since 
degree R(x) < m — 1, using Lemma [2] we get 



R{x)V-\x) = Y,— 



AC— 1 

where the matrices C m ± are given by 

C ^ k = <a TrTr 1 sW? ^OwX-Mi (V r ro ( a 0))< , *- 1 >( a ; m , fc ). 

(det {Vm{x))) (lk > {x m , k ) 

According to Lemma [2 V m {x m ^) (Adj (KhOb)))^* -1 ^ (»m,fc) = OjvxAf and taking 
into account that R(x m> k) = V(x m ^) — C{x m ^)V m {x m ^), the previous expression 
becomes 

C ^ k = (a *<v ( ^W? T^ m , fe )(Adj (^(x)))^- 1 '^^)- 

(det (K m (a:)))^J(:r m ,fc) 

Then, 

V m {x) 



V(x) = C{x)V m {x)+Y,C; 



m,k ' 

•E Xm.k 
k—1 

Since 
we have 



V(x) = C(x)V m (x) + ^2c 
Taking x — h(t), we have 



rn , k " 

•Em.k 

k—1 



v(h(t)) = c(h(t))v m (h(t))+i2 c ^ Vm{Xm,k) lz\ m - 

X mik ~ h[t) 

Then, from the integral representation of the associated polynomials of the first 
kind 

B ™- i(z)= y — v^htx) — dw ( h w> 

it follows that 

V(h(t))dW(h(t)) = J C(h(t))VUm)dW(h(t))^^C m ^l 1 (x m ^ k ). 

k—1 

So, from the orthogonality of {Vm}meN with respect to W we have 

f V(h(t))dW(h(t)) = ]T C^.fcSS-l^m.fc), 
^ k=l 

and the statement follows. □ 
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The next result is an extension of one proved by A. J. Duran in [5]. It deals with 
the ratio asymptoticss of the m-th orthogonal polynomial V m with respect to the 
generalized Markov matrix function, 3", and the (to — l)-th associated polynomial 

of the first kind T>^_ v 

Theorem 16 (Generalized Markov's theorem). Let 11 be a quasi- definite vector of 
linear functionals, {V m } m £f$ be the sequence of matrix polynomials left- orthogonal 
with respect to the generalized Markov matrix function, 3, defined by (|22|) , and let 
\fS>m } TO eN be the sequence of associated polynomials of the first kind for {H>m}meN 
and VL. Then, 

lim V m 1 (z)^_ 1 (z) = n.z) 

m — >oc 

locally uniformly in <C\T, where T = (~]n>qMn, Mjv = U„>jv{zeros of V m } . 



Proof. First, from Lemma [5] we get 



v- 1 (z)v£l 1 (z) = J2r m>h 



k=l 



%m,k 



where r m .fc are the matrix coefficients that appear in the quadrature formula pre- 
sented in Theorem 1151 and x m .k are the zeros of V m . On the other hand, there 
always exist complex numbers y m ,k such that /i(j/ m ,fc) = Xm,k, and 



i 



fc=i 



z - h(y mjk ) 



We consider the sequence of discrete matrices of measures {/i m } m gN defined by 

Mr, 

Thus, 



k=l 



i 



' z - h(y mik ) 



dfi m (h(x)) 
z — h{x) 



(34) 



if z is not a zero of V m . Taking into account (|34[) . it will be enough to prove that 

dfim(h(x)) 



lim 



oo J z — h(x) 



= J(z) for zec\r. 



The first step deals with the pointwise convergence. Otherwise, we assume that 
there exists a complex number z G C \ T, an increasing sequence of nonnegative 
integers {mi}i^n, and a positive constant C such that 



dn mi {h(x)) 



■J{z) 



> C > 0, I > 0, 



(35) 



z — h(x) 

where || . ||2 denotes the spectral norm of a matrix, i.e., 

\\A\\ 2 = max {VA: A is a eigenvalue of^4*v4}. 

Taking an increasing sequence {a k }keN such that a k — > oo, and using the Banach- 
Alaoglu's theorem there exists a subsequence {n}; e N from {mi}i^, defined on a 



MATRIX ORTHOGONALITY ON THE REAL LINE 



25 



curve 7^ contained in a disc \z\ < at, with the same fc-th moment of the vector of 
linear functionals, U, for k < 1r\ — 1, such that 



lim 

/ — >oo 



f(h(x))dfx n (h(x)) 



1 



y (^) 

; — h(x) 



dz . 



(36) 



Moreover, 

f dnn(h(x)) 
J z-h(x) 



J(z) 



< 



Ik 



dfi ri (h(x)) 
- h{x) 



djjL ri {h(x)) 



h{x) 



with £k in the exterior of the disc \z\ < at- We write So for the first moment of the 
matrices of measures \i Tl which is the first moment of U. Then, by taking k and 
then ri large enough, from (|35p and (|36p we obtain 



C 
2 



< max 



< max 



1 



h(a k )\ 
1 



dn n {h{x)) 



IS', 



2 • 



\z - h(a k )\, 

But this yields C = and, therefore, ([3"5]) is not possible. 

The next step is to prove that the analytic functions which are the entries of the 

matrix [ < ^^ m ^ Y 1) are uniformly bounded in compact sets of C \ T. Then, the 
J z- h(x) 

uniform convergence in compact subsets of C \ T will follow from Stieltjes-Vitali's 
theorem. 

Given a compact K C C \ T, let notice that K fl Mm ^ 0, for N large enough, 
and then there exists A > such that 



1 



h{x) 



Then, for n > N 



The spectral norm 



< A, for z e K and h(x) e Afjv. 
d/i n (/i(a;)) 



/ 



<^Sn 



d/j, m (h(x)) 



z — h(x) 

equivalence of the norms in finite dimensional spaces, the result follows. 



is uniformly bounded and, therefore, from the 

□ 



Remark . In an analog way we can deduce the following result. Let {G m } me N be 
the sequence of matrix polynomials right-orthogonal with respect to the generalized 
Markov function 5F and let {Gm } m gN be the sequence of associated polynomials of 
the first kind for {G m } m( =N and U. Then, 

lim G%_ 1 (z)G- 1 (z)= ?(z), 

m — >oc 

for zeC\T and the convergence is locally uniformly on C\T, where 
T = n N>0 M N , M N =~ T 



>N 



{zeros of G m }. 
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